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Extra-large deformations in ultra-soft elastic materials are ubiquitous, yet systematic studies and methods to 
understand the mechanics of such huge strains are lacking. Here we investigate this complex problem systemat¬ 
ically with a simple experiment: hy introducing a heavy head of radius a in an incompressible ultra-soft elastic 
medium. We find a scaling law for the penetration depth (S) of the bead inside the softest gels as <5 ~ 
which is vindicated by an original asymptotic analytic model developed in this article. This model demonstrates 
that the observed relationship is precisely at the demarcating boundary of what would be required for the field 
variables to either diverge or converge. This correspondence between a unique mathematical prediction and the 
experimental observation ushers in new insights into the behavior of the deformations of strongly non-linear 
materials. 

PACS numbers: 46.25.-y,83.85.Cg,46.05.+b,83.80.Va 


Singularities are pervasive in various problems of linear continuum mechanics. In wetting, stress diverges at a moving contact 
line mill; it diverges at tbe tip of a crack or even at a sharp point indenting a plane fJl. Understanding how such singularities 
can be tempered has often given rise to new physics invariably prompting us, on many occasions, to investigate a material 
phenomenon at a molecular dimension and then herald a way to bridge the near field with the far field behavior in a rather 
non-trivial manner. Nature, however, performs the difficult task herself and leaves her signature in a way that is independent 
of the constitutive property of a material, yet it belongs to a class of universality. What we report here is such a universality 
that is discovered in the large deformation behavior of ultra-soft gels. Soft solids undergoing huge deformations exhibit various 
fascinating and puzzling mechanical behaviors iHISl. Our experimental protocol to study extra-large elastic deformations is 
remarkably simple, in that a heavy bead of stainless steel is gently deposited on the horizontal flat surface of a gel. The compliant 
gel is deformed by the load exerted by the heavy bead. It reaches a stable (elasto-buoyant) equilibrium position when the elastic 
force exerted by the sutTounding gel balances its weight m, within few tenths of a second. This experiment can be viewed 
as an elastic analog of the falling ball viscometry, in which the bead reaches a terminal sedimentation velocity resulting from 
the balance of the bead’s weight and the viscous drag force M- In the limit of Hookean elasticity, an analogy with the Stokes 
equation (by replacing shear viscosity with shear modulus and velocity with depth of submersion, suggests that S ^ a^, 

a being the sphere radius. While for the higher modulus gels (Figj^a), such a relationship is more or less valid, for the softer 
ones, when the bead is totally engulfed by the gel and the deformations are very large (Fig{T]-b), it is observed that the depth 
scales with the bead’s radius with an exponent of 3/2. This is a non-trivial result that cannot be explained by the usual neo- 
Hookean model, i.e, by considering a quadratic elastic energy density with respect to finite strains. According to this model, the 
displacement would be infinite! When the same problem is analyzed using a new analysis presented in this paper, it is found, 
remarkably, that the exponent of the radius of the bead is just at the juncture of what would be required for the field variables to 
avoid divergence. What is quite remarkable about this analysis is that no detailed non-linear behavior of elastic solids MM 
needs to be speculated. The scaling laws are universal that being independent of the particular constitutive law of the elastic 
material. This unique correspondence between a mathematical prediction and the experimental results, thereby unfolding new 
physics of highly non-linear deformations, is the subject of this paper. 


I. EXPEIUMENTAL EVIDENCE OE ELASTO-BUOYANCY 


Cross-linked polyacrylamide gels are used in the experiments reported below. The gel solutions are prepared according to 
the recipe published previously dEa and cured in home-built glass containers (70 mm x 50 mm x 40 mm). The inner walls 


serge .mora @ umontpellier.fr 



2 



FIG. 1: Snapshots and schematics of the experiment. (Top) Side views of two transparent cells filled with a polyacrylamide gel having a shear 
modulus of 1160 Pa (a) and 13 Pa (b). Two identical steel beads (5 mm diameter) have been deposited in the free air-gel interface. The vertical 
downshifts are respectively S — 0.03 5o and 5 = 320 Sq. (Bottom) Schematic of the experiments. 


of the containers are grafted with a thin layer (~5 nm) of polydimethyl siloxane chains so that the gel solution contacted the 
walls at 90° to ensure that the surface of the cured gel is flat. All the experiments are performed after two hours of gelation. For 
the estimation of the shear modulus of the gel, we use a linear elastic model in order to ensure consistency of quantification for 
all the gels. The shear modulus is determined from the resonant mode of vibration of a gel slab confined between two parallel 
glass slides Ea. Steel spheres (density 7.8 g/cc, diameters 1-10 mm) are gently placed one by one on the gel surface and its 
side-view image is captured by a camera. Dissipative processes within the gel dampen any oscillations, and the spheres sink until 
they become stagnant in the polyacrylamide gel, the whole process taking only fractions of a second. The depth of submersion, 
S, is measured from the upper surface of the gel till the base of sphere, that denoting the net downward displacement due to 
the inclusion of the spherical particle by the surface (Fig{T]l. The cells are large enough to avoid any finite boundary effects, 
including the side walls and the bottom. The measurement of depth for each sphere is made in the central region of the container. 
After each measurement, the sphere is gently removed from the gel using a magnet, held slightly away from the free surface. We 
wait for a few minutes between each measurement that ensures there is no memory of the position of the previous sphere inside 
the gel. 

If the bead is too small and the gel is stiff ||7] [^, the surface bends slightly under the weight of the bead and S a 
(Figj^a). By increasing the bead radius or decreasing the elastic modulus, the particle submerges itself to a considerable 
depth inside the gel. The surface of the gel wraps around the particle and closes to create a line singularity connecting the 
particle to the free surface of the gel (Fig[T]-b). Strings of tiny air bubbles appear in this thin channel, which soon coalesce 
and escape through it while the channel further closes due to the auto-wetting forces of the gel’s surface. If the surface of 
the gel is premarked with ink spots, it is easy to visualize that the surface of the gel becomes appreciably stretched while the 
sphere sinks through the gel while being still connected to the free surface via a thin channel. It is also possible to release ink 
inside the gel in the form of thin vertical lines with the help of a fine needle, which bend toward the sphere in a dramatic way 
when the sphere is released inside the gel demonstrating a substantial amount of tensile strain developed in the gel network 
due to the inclusion of the bead. These basic experiments were reported in a previous article M, but without a detailed anal¬ 
ysis. Here we report a detailed set of experiments, in which the shear modulus of the gel was varied from 13 Pa to about 3000 Pa. 

Prior to subjecting these experimental results to a comparison with a theoretical analysis, we needed to verify how meaningful 
it is to consider only the effect of elasticity by ignoring the surface tension of the gel in predicting the depth of submersion of 
the bead and how reversible is the deformation of the gel. First question is partly philosophical that rests upon the distinction 
between surface free energy 7 and surface tension. The latter differs from the former by surface stress dj/de, e being the 
surface strain. As the major constituent of these amorphous gels is water, we expect that the surface stress is negligible. In fact, 
several recent studies that measured surface tensions of various amorphous soft polymers strongly suggest that their surface 
tensions are practically same as their surface free energies EM . Thus we need to figure out only if the surface free energies 
of the gels play any role. Part of the answer can be obtained by comparing them with the equivalent spring constant of the 
sample. The latter can be obtained by slightly raising the height of the bead by an electromagnet and releasing it so that the bead 
undergoes an under-damped oscillation and reaches the neutral position. For three gels, in which the beads were completely 
submerged, the equivalent spring constants were estimated from the frequency of oscillation to be 0.2 N/m, 5.6 N/m and 13 
N/m, which increase systematically with their shear moduli (13 Pa, 140 Pa and 360 Pa). Fig|^a shows a typical profile of such 
an oscillation. Comparing these spring constants with the surface tension (0.07 N/m) of water, we conclude that the contribution 
of surface tension can be safely neglected for all the gels used in this study except, perhaps, for the lowest modulus gel for 
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which the spring constant is three times that of the gel’s surface tension. However, when a bead is completely submerged in the 
gel, any variation of the height of the sphere does not alter the area and the excess energy of the free surface of the gel. We 
thus believe that the surface tension spring does not play a significant role in determining the depth of submersion of the sphere 
as long as it is completely engulfed by the gel. Further support to this viewpoint, i.e. the dominant role of the elasticity over 
surface tension was gathered from the experiment described below that also exemplified the reversibility of the softest gel (p ^ 
13 Pa) employed in our experiments. 

By varying only the system temperature, the elastic modulus of polyacrylamide gels changes, which increases with tempera¬ 
ture. A thin layer of Paraffin oil is poured over the gel to avoid its surface from drying. The temperature of the gel is monitored 
by placing a thermometer inside an identical sample of gel in a similar sized container, placed inside the oven. After the gel 
is heated to 70 °C, a 5 mm diameter steel sphere is released into it through the layer of paraffin oil. The depth of the sphere is 
measured at this temperature while it is in the oven. As the gel is gradually cooled, the sphere sinks deeper inside it. We wait 
for an hour between acquiring data for the depth of sphere at each temperature to allow the gel to equilibrate reasonably well. 
After the gel is cooled to about 5 °C, it is heated again that decreases the depth of the embedded steel sphere (Heating Cycle, 
Fig|^b). The depth of submersion of the sphere plotted as a function of the temperature for both the cooling and the heating 
cycles (Fig|^b) shows that there is a little hysteresis in this system, in that the difference in the depths of the sphere for a given 
temperature is within 5%. We conclude that the deformations of the gel generated by a bead are predominantly reversible and 
controlled by elastic forces. 



t (s) Temperature (“C) 


FIG. 2: Experimental evidence of reversibility of the gels, (a) A 10 mm diameter steel sphere immersed in 140 Pa gel, when slightly disturbed 
from its elastobuoyant position via an electromagnet, undergoes under-damped oscillations about its equilibrium depth 5eq- (b) Depth of 
submersion 5 of a 5 mm diameter steel sphere in a soft gel (p ~ 13 Pa) varying as a function of its temperature. The experiments in the cooling 
cycle were performed first following which the gel was heated systematically to obtain the data for the heating cycle. 


The penetration depth, 6, is plotted versus the radius a of the steel spheres in log-scales in upper inset of Figj^ for ten 
different gels (/r ; 13 Pa - 2930 Pa) and different radii. The upper white region of this inset shows the data points for beads that 
were completely below the surface of the gel. When we non-dimensionalize the depths as well as the radii of each bead in each 
gel by the material lengthscale Jq defined as Jq = where fi is the shear modulus of the gel, Ap is the effective 

density of the buoyant spheres, and g is the gravitational acmeration, we see that all the data cluster around a mean master 
curve with two distinctive asymptotic limits. For each sphere-gel system, Ap was precisely estimated from the image on the 
basis of the Archimedes principle of buoyancy analogous to the scenario for floatation on a liquid ll25l . The normalized data 
can be divided more or less into two regimes; one with the non-dimensional radii q/Sq < 1 and other one with a/So > 1 with 
an intermediate transition regime. We fit the non-dimensional depths as a function of the non-dimensional radii for the regime 
a/6o < 1 with the power law function {6/6q) = k{a/6o)°' with adjustable parameters, a and k and find a = 1.96 ± 0.06 and 
k = 1.09 ± 0.1. The error bars are obtained from a 95% confidence limit analysis. We conclude that for a/6o < 1, the depths as 
a function of the radii follows 5 ^ a? within the error limits. On careful examination of the experimental points corresponding 
to large deformations, i.e. a > Sq, we find that the data corresponding to the two softest gels (p : 13 Pa and 25 Pa) are shifted 
from the rest of the data due to their multiplicative factor {k) being significantly larger than the rest of the data for the other gels. 
This indicates that the lengthscale Sq, which is defined with the elastic linear properties of the sample, is therefore not sufficient 
to describe the whole data accurately, and non linear effects have to be taken into account. We conclude that all the data in 
this regime (a/So > 1) cannot be investigated together. It is more appropriate to investigate the data for each gel composition 
separately i.e, for a given non-linear material stress-strain relationship. By fitting the data for the two softest gels {p, : 13 Pa and 
25 Pa) where all the beads are completely engulfed, with the power law function, we find that ai 3 pa = 1-42 ± 0.05 (see lower 
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inset of Figj^, fci 3 Pa = 2.52 ± 0.3 and 025 pa = 1-52 ± 0.09, ^25 Pa = 1-63 ± 0.3. The exponents for the fits for the gels 
(a/(5o ~ 1) in the intermediate regime lie between a ^ 1.5 and a ^ 2. Thus, from the experimental observations, we infer that 
in the limit of a significantly greater than 60 , the general trend is close to S ^ a^, where p is in the range of 1.4-1.5, within the 
error limits. Thus, we conclude that the power law observed (<5 ~ a^) in the gels of higher shear moduli is closer to the regime 
already studied before what one would expect with an analogy with the elastic Stokes equation in the Hookean limit. What is 
astounding is the observation of J ^ in the case of the gels where S > 2a for each tested bead radius. In order to interpret 
these novel observations in the ultra-soft gels, we develop a new model to tackle such extra-large deformations in the following 
part of the paper. 
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FIG. 3: Dimensionless depth of spheres (S/Sq) plotted as a function of its dimensionless radius a/So, for various shear modulus of the gels 
from /r = 13 Pa to 2930 Pa. The grey lines indicate power law curves (S/So ( 0 /^ 0 )“) with q = 2 and a = 1.5 for the two asymptotic limits 
for the normalized data. (Upper Inset). Depths (ti) versus radii (a) for all the spheres. The plot area is divided into two domains; the boundary 
indicating S = 2a. The data points above the boundary indicates that the spheres were entirely below the gel’s surface. (Lower Inset). Best fit 
(5/5o = k{a/So)°‘) for the softest gel (13 Pa) highlighted (see text). 


II. ASYMPTOTIC ANALYTIC MODEL FOR LARGE ELASTIC DEFORMATIONS 


From the outset, a motivating picture of the problem can be gleaned from the comparison of the potential energy of the bead 
in the gravitational field and the energy of the elastic deformation of the gel. In the limit of Hookean deformation, an analogy 
with the Stokes equation suggests that the elastic energy is of the form pS^a, which is to be compared with the gravitational 
potential energy Apga^S. This leads to the scaling: <5 ^ a^. However, the experimentally observed scaling S ^ implies 
that the elastic energy in extremely large deformation must scale as 6 ^. The detailed analysis based on a model, presented in this 
paper for the first time, shows that the above scaling is non-trivial. Moreover, the scaling result that follows is independent of 
the constitutive laws of an elastic material. 


A. Gravity energy of engulfed spheres 

Consider deformations that preserve the volume, as is the case with the elastic gels used in the experiments. Since the bead 
is supposed to be totally engulfed with 6 a, the surface of the deformed gel is fairly flat and horizontal. A supplementary 
(virtual) downshift 6 ' produces an opposite rise of the same volume of gel. Therefore the gravitational energy variation of the 
system consisting of the sphere plus the gel is ^na^ApgS'. We conclude that the gravitational energy shift is £gr ~ ^ira^ApgS 
in the limit we are considering. 


B. Outline of calculating the elastic energy In the limit 5 ^ a 

The axis of symmetry being the vertical axis, the two coordinates changed by the deformation are the radius in the horizontal 
plane, r, and the vertical coordinate, z. The deformation maps the undisturbed state with coordinates (r, z) to a disturbed state 
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(i?(r, z), Z{r, z)), or, in radial coordinates from coordinates (f; 9) to {R{r, 9); Z{r, 9)) with f radius and 9 polar angle of the 
(r; z) plane (see Fig|^(d)). The elastic energy is a function of the strain tensor (also called right Cauchy-Green deformation 
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FIG. 4: Sketches corresponding to the steps of the calculation of the elastic energy, (a) Reference state (no deformation); (b) A point-load is 
applied at the free surface. The displacement at the application point is <5. (c) A sphere of radius a indents the free surface over the distance 
S. (d) Mapping from the reference state (solid lines, with the bead at the surface) to the deformed state. A point of the gel in the rest state is 
located with coordinates (r, z). f is the distance from the initial contact point of the bead. In the deformed state, the point that was at (r, z) is 
located at R{r, z), Z{r, z). 

tensor), which gives the square of local change in distances due to deformation: C = F'^F, F being the deformation gradient 
tensor ll26l . In the absence of any preferred direction, the elastic energy of an incompressible solid may depend only on two 
scalars (invariant under global rotation) that can be made out of the strain tensor: the trace (Ii) and the sum of the square of its 
components (/ 2 ). In cylindrical coordinates with an azimuthal invariance, these invariants read Ii — R\ + ^ + R^^ + Z\ + Z\ 

and I 2 = {R\ + Z\Y + + ‘^{R,rR,z + where indices preceded by a comma denote partial 

derivatives. For the neo-Hookean model, which is commonly used to describe soft gels, W = ^/i. However, beyond a certain 
deformation, the neo-Hookean model cannot be a fair representation of rubber-like materials whose elasticity originates from 
unfolding of polymer chains: once these chains reach their full extension, the energy cost for a supplementary unfolding 
diverges so that a further deformation is accompanied by a diverging additional elastic energy, in contrast with the ideal 
neo-Hookean law (see Fig. |^. In what follows, we do not restrict ourselves to the neo-Hookean case. 

The deformation field of the gel can be described with two characteristic lengths, the downshift of the sphere 8 , and the radius 
of the sphere a. Here, we are dealing with the limit of large downshifts, that is 8 far larger than a. Below we demonstrate that 
in this asymptotic case, the elastic energy of the system does not depend on a in the limit a 8. As a first step we consider 
a normal point-load applied on the free elastic surface at (r, z) = (0, 0) (Fig|^(b)) and we find a scaling law for the elastic 
energy per unit volume at a given distance f from the loading point. The condition for the elastic energy to be convergent is 
established. Then as a second step we consider a hard sphere with a finite radius, producing a displacement 8 at the contact point 
(r, z) = (0, 0) equal to the displacement at (r, z) = (0, 0) induced by the previously considered point-like load. We demonstrate 
that for (5 <C a, the elastic energy associated with the displacement field generated by the sphere or the point-like load are equal: 
These are independent of a and proportional to 8 ^. 
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C. Elastic energy for a point-like load 

Incompressibility of the elastic medium is imposed by writing that the determinant T) of the first derivatives of R(a;, y, z) is 
equal to one. For axisymmetric deformations, V = ^ — R,zZ^r)- The total energy of the system {gel + sphere} reads : 


4 pO POO 

£ = -TrApa^gZ{0,0) + 2Tr / dz dr{W — qT>)r, 
3 J-oo Jo 


( 1 ) 


where the Lagrange multiplier q{r, z) imposes the incompressibility condition T> = 1 lIZTl . 

For the point-like heavy sphere we first consider (see Fig. |^(b)), the weight the bead can be reduced to a point force located 
at (r, z) = (0,0). Eq. [^simplihes into: 


£ = 


IGtt^ 


/ O poo pO poo 

dz / r2'(r)i5^(r)dr-I-27r / dz dr{W — q'D)r, 

-oo Jo J — oo J 0 


( 2 ) 


where 6^{r, z) is the 3D Dirac distribution. 

The Euler-Lagrange conditions of minimization of the energy (Eq. read 0 : 


d (W — qV) r 
dR^r 

5 (W - qV) r 
dZ~r 


+ 


+ 


a (W - qV) r 

d{W- qV) r 

Wz 


d (W — qV) r 
87rrApga3 ^3^^^ 


(3) 


(4) 


Let 5 be the displacement of the gel at (r; z) = (0; 0). ^ is the unique relevant length for the displacement field. Eor points located 
at a distance to the load f much smaller than <5, S is no more relevant at these short length-scales. Therefore no length-scales are 
expected to occur in the scaling laws and one assumes a power law of f for the displacements for r 6: Z = 6 + r^fi{d) and 
R = f 2 ( 9 ). Since the gel is vertically stretched and horizontally squeezed in the vicinity of the bead one assumes that 7 > /3. 
This hypothesis will be checked afterwards (see Eq. [T 6 ]below). 

The constraint of conservation of volume reads (for f <C 5): 




R 

r 


2? = 1 => — ( R J, 

The condition that 'y > (3 with Eq. [^yields /?< 1 . 


Z,z — R,z Z^r ) = 1 ^ 7 = 


3-/3 


(5) 


We assume that the strain energy density function is W ^ (the exponents ai and a 2 are constant) in the range of 

strains undergone by the elastic solid in the vicinity of the application point, i.e. for f far smaller than S. This choice for W is 
crucial to obtain the scaling law, but it has no effect on the hnal results, as argued in section IIE it does not limit the general 
nature of the theory. 


The scaling laws for the hrst and the second invariants of the Cauchy deformation tensor are: 


R^ 


f.2fS-2 

7 I 


I\ — Tr{C) — R ^ ^ R z: Z z- + Z 


,f2P-2 


f:2/3-2 


r 

72 ^2 




J2 = Tr(C^) = (i?", + Zt) + {Rz + Z%Y + 2{R^rR,z + Z^zZ,rY + — 


3-4 


(6) 

(7) 


These scaling laws for 1\ and I 2 are used to obtain the scaling laws for the various terms of the Cauchy-Poisson equations (Eqs. 
|^and|^. We hrst deal with Eq. 


diW- qV) r dW dh dW dh dV 

= -b J'TT.—4:714- qT -, 


dR , 


dll dR^r dl2 dR^r 


dR, 


(RriRl + zl) + R,z{R,rR,z + z^zZA -q RZ,z 
oil 0 I 2 ’ ’ 


^r^+2(c«i+2c,2-1)0-1) 


.,f;^+2(c.i+2c2-l)(/3-l) 


3+1 

2 


( 8 ) 
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In the same way one finds: 


d{yV-qV)r 3-/3 


+ 2(ai + 2a2-l)(/3-l), 


(9) 


d{yV — q'D)r .. g-i 
dR 


( 10 ) 


From Eq. |^and Eqs. [^|^and 
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one obtains f 2 ^+ 2 (“i+ 2 a 2 i)iP i) ^ gf ^^2 , and the scaling law for the Lagrange multiplier: 

(/3—l)(2ai+4a2 —3) (11) 


q ^ r 

The last Cauchy-Poisson equation (Eq. is now used to get an expression of /3 as a function of ai and a 2 . It writes: 

,3 


I f d (W — qT>) \ ^ f d (W — qT>) ^ _ SnApga ^3 


- I r 
r 


dZ, 

djW-qV) 


dZ, 


5\v). 


( 12 ) 


Denoting as azr = ~ ^^''dz quantities present on the left-hand side of Eqj^ one finds that this 

left-hand side is the divergence (expressed in cylindrical coordinates) of the z-components of the stress tensor a 
We first calculate the scaling laws for a^r and a^z using Eqs. |^|6 11 


a (W - qV) dW dh dW dh 


dZ, 


dh dZr dh dZ, 


dV 


dyv R 

2 —Z, +q-R,z. 

dll r 


(13) 

(14) 


..,f(/3-l)(2<»l+4Q2-l) 


pl-3 


The second term of the right-hand side of Eq. [^scales as i)( 2 Qi-i- 4 a 2 3 )^i ^ p negligible when f tends to zero when 
compared to the first term of Eq. provided that the exponent of the second term is larger than the exponent of the first one, 


i.e.: 


(/3 - l)(2ai + 4^2 - 3) + 1 - /3 > (/3 - l)(2ai + 4^2 - 1), 


(15) 


which is formally equivalent to /3 < 1. Since it is assumed from the beginning that /3 < 1, we conclude that Uzr 
f(/ 3 -i)( 2 ai+ 4 a 2 -i) jjj jjjg g^nie Way, (Tzz ~ ^(/ 3 -i)( 2 ai+ 4 a 2 -i). ^follow the Same power law when f tends 
to zero. The right-hand side of Eq. f^being a delta-like charge density, we conclude from Gauss’s theorem that the scaling law 
for Gzr and Gzz is also azr (Zzz Therefore, f(/ 3 -i)( 2 ai-i- 4 a 2 -i) ^ ^- 2 ^ 


2ai + 4a2 — 3 
2ai + 4a2 — 1 


(16) 


/3 is an increasing function of (ai + 20 : 2 ). It is negative for ai + 2a2 < 3/2, and it is always smaller than 1. The displacement 
at the application point is finite if the exponents are positive, i.e. if ai + 2 a 2 > 3/2. Otherwise the material cannot withstand a 
point-load. In the following we assume that ai -f 2 a 2 > 3/2, an ansatz that will be legitimized later. 


Since 8 is the unique length-scale of the deformation, the coordinates in the deformed configuration can be written as: R = 
i5/o (I) and Z = Jpo (!)> where /o and po ure two numerical functions depending only on the constitutive law of the elastic 
medium. The invariants Ii and h, and thus the elastic energy density are dimensionless functions depending only on r/S. The 
elastic energy due to the deformation, 

£eio= fff Wod^r, (17) 

JJJ 2 > 0 ;f >0 

(subscripts 0 are for the point-load problem) is therefore proportional to 6^ with a coefficient Cq, which has the dimension of 
shear modulus and depends on the mechanical properties of the elastic solid, as: 

Eel 0 = CuS^- (18) 

Note that the convergence of the integral in Eqj^is ensured since the work done by the applied force is finite for ai + 2a2 > |, 
yielding the validity condition of Eq. [T^ 
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D. Elastic energy for a finite sphere 


We consider now a heavy bead of radius a, and denote 5 its vertical downshift, which is supposed to be far larger than a (Fig. 

§(c)). 

• The unique relevant length-scale for f 3 > a being S, the displacement field reduces in this range of f to the disp lacement 


field for a point-load with the same penetration depth 6 : Z = SgQ where po has been defined in section 


lie 


• The power law for Z — S with the exponent given by Eq. [T^ applies in the intermediate range a f S. Combining it 
with the previous expression of Z (valid for f ^ a) yields: 


Z-6r~.5 




• Close to the bead if ^ (5) the unique relevant length-scale is a. Z — S can be expressed as 

'r' 


Z-S = Kg, 


(19) 


( 20 ) 


where gi is a numerical function and the constant K (with the dimension of a length) depends on a and also on the far field 
deformation, i.e. on 6 . Comparing Eq. 19 with the general expressions for Z in the range f 6, one obtains K = 6 (f )^- 

A similar expression for R can be obtained with the exponent (3 — /3)/2, leading to negligible contributions for the elastic 
energy. 

Eor f S, the first invariant scales as Z\ and Z\, and the second invariant l-i scales as Z\ and Z\. The strain energy 

a\i 

sJ 
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with the expression of K) as W 


density function W ~ can thus be written for f S (using Eq. 

In the range f 3> a, W — Wq, Wq being the elastic energy density of the point-load problem (Eig. 1^. 

, . R _ 2 / \p _ 3 

Moreover, in the range f one finds from the scaling laws of section '''' 


3-3 




lie 




3-3 



FIG. 5: Sketch of the energy density profiles Wo (e.g. dotted line) and W (e.g. solid line) as a function of the distance t f to the initial contact 
point. The dashed horizontal double arrows highlight two domains, corresponding to (i) f ^ a where the elastic energy densities W and Wo 
are similar, and (ii) f <C 5 where they are different. 


One concludes that, for any f, W — >Vo (f)^ ^53 (q)i with g^ a numerical function (depending neither on 6 nor on a) 
whose limit as its argument approaches oo equals 0. Thus, we obtain: 

/r»\ _/T*\ / n \ d 

( 21 ) 




18 


that the elastic energy 


Since one assumes again that ai -F 2 a2 > 5 , /3 > 0 and the difference of the elastic energy with the bead of radius a, £ei, to the 
elastic energy with the point-load, 0 , is negligible with respect to 5^. One concludes from Eq. 
is proportional to S^. 

Note that if ai -F 2 q ;2 were smaller than 3/2 the near-held part of the elastic energy would not be negligible anymore compared 
to the far field contribution. 
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E. The strain energy density function 


In sections 


lie 


and 


IID 


the strain energy density function has been assumed to scale as W ~ for the strains 

encountered at f ^ d. It has been then demonstrated that the elastic energy corresponding to the downshift of a sphere of finite 
radius a over the distance 6 is proportional to 6 ^ in the limit J 3 > a if ai + 2 a2 > 3/2. 


It is worth considering the case of an elastic material for which the strain energy density function increases softer than 

with ai + 2 a 2 = 3/2. For f S we demonstrated that Z — S = S (^) and then ^ showing that the 

strain is arbitrarily large in the limit S ^ a , if /3 < 1, i.e. if ai + 2a2 < i/2. However, due to the finite maximum stretch 
of the polymer chains constituting the material, the maximum stretching of any real elastic rubber-like material is bounded. To 
this maximum stretching corresponds a divergence of the strain energy function. This divergence results in a steeper and steeper 
increase of W with Ii, I 2 , or both, which is associated with increasing values of the exponents ai and/or a 2 lf26l . 

In order to illustrate this last point, let us consider a particular strain energy density function, for instance the Gent 
hyper-elastic model ETl l28l . The strain energy density function of this model has a singularity when the first invariant Ii 
reaches a limiting value. It is plotted in Fig. deviations from the initial neo-Hookean behaviour (ai = 1 and a 2 = 0) 
yield an increasingly stiffer and stiffer strain energy density function. At any value Iiq of Ii one can define the local exponent 
0^10 = ^ so that at the vicinity of Iiq the strain energy density function is W ^ locally. Beyond a certain finite value 
//q of /lo the exponent is larger that 3/2, ensuring that for any Ii larger than J/g the condition ai + 2a2 > 3/2 is fulfilled. 
Note that values of Ii larger than //g are automatically reached, otherwise the strain would diverge around f ^ 0, leading to 
arbitrarily high values of Ii (as explained above). 


What is illustrated using the example of the Gent model is general, and can be applied to any elastic constitutive law involving 
Ii and I 2 of a real elastic material: starting from low to moderate strains for which the neo-Hookean model is expected to apply 
far from the bead, energy density functions stiffer than with ai + 2 a 2 > | are necessarily encountered next to the bead 

for any real elastic material. 

A material for which the failure limit is reached before the strain energy density is stiffer than with ai + 2a2 > i/2 

would not be able to sustain the heavy sphere, and thus would be drilled by the sphere. These cases can be definitively excluded 
for the experiments we are dealing with since no fracture, plasticity nor creep have been observed during these experiments. 



h 


FIG. 6: Strain energy density function given by the incompressible Gent material model W = log ~ ^ j im with Jm = 97, 

plotted in log-log scales (black solid line). ,Jm -F 3 is the limiting value of the first invariant I\. The strain energy density function for the 
neo-Hookean material is plotted with dashed line. The two gray straight lines indicate slopes 1 and 2, i.e. values of the local exponent ai equal 
to 1 and 2. 7 *q is the value of the first invariant beyond which the local exponent is larger than 3/2. 

One concludes that for any elastic material, the elastic energy of the gel is proportional to 5^ within the limit 5 ^ a. The 
proportionality constant is Co (see Eq. [T8]l depends on the whole elastic constitutive law of the elastic material. 
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F. Equilibrium condition 

The energy is minimum in the equilibrium state. The i5-derivatives of 8 gr and 8 ei are therefore equal giving the scaling law 
valid for S ^ a, for a given elastic material; 


5 ~ 0^/2 {Apg)^/^ . (22) 

Note that the beads has been assumed to be spherical although this is not a crucial point in this theory: if it is compact with all 
typical lengths of the same order a the previous scaling arguments apply. 


III. DISCUSSION AND CONCLUSION 

The theoretical analysis presented above provides an understanding sui generis of the experimentally observed scaling law. 

The elastic material being chosen (so that the proportionality constant Cq is fixed), this law predicts the scaling behaviour 
of the downshift of spheres with different radii and/or different densities. These predictions are in quantitative agreement with 
the observations; one observes experimentally a scaling law with an exponent close to 3/2 for the radius, and the prefactor 
is not only related to the linear elastic properties (/r), but also to non-linear features of the elastic material. These non-linear 
features being a priori distinct from one gel to another one, this explains why plotting 5/5q as function of a/5 q for different gel 
compositions yields to different curves. The prefactor in the scaling law depends not only on the shear modulus of the solid, 
but also on the elastic behaviour at large deformations, in agreement with the experiments (HD). This provides a way to assess 
some characteristics of elastic materials under large strains, as strain stiffening properties. 

The range of applications of the derived scaling goes far beyond the description of the elasto-buoyancy phenomenon. For 
instance, it can be directly checked that the recent indentation experiments performed on various types of compliant gels, by 
Fakhouri et al ua, are in quantitative agreement with our prediction, however their very naive model based on the neo-Hookean 
equation, misses the driving non-linearities. 

In this paper, we have shown that the exponent 3/2 for the radius is independent of the strain-stress relation provided that 
the increase in the elastic energy density (W) with the strain is stiff enough. Generic behaviors for elastic materials undergoing 
large and complex deformations can therefore be identified, going beyond scaling arguments that is blind to the crucial effect 
of strain stiffening. This work sheds new light on the mechanics of extremely large elastic deformations that are crucial in 
various emerging techniques that even encompass an important procedure such as the computer-assisted surgery involving human 
organs IIIHSII, which are indeed non-linear materials and undergo large deformations. In all these important fields, an ideal 
neo-Hookean model is used that, according to our current work, is not valid sensu stricto. The danger is that one may miss 
the divergences of the solution that is inherent in the ideal neo-Hookean model because the full solution is not used. The 
work presented here is a panacea to the possible pitfalls that one may embark upon in studying extreme deformations of non¬ 
linear materials, thereby also offering possible benchmarks for numerical simulations and even opportunities to formulate new 
simulations methods. 
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